Pseudogap Value in the Energy Spectrum of LaOFeAs: Fixed Spin Moment Treatment 
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The experimental data available up to date in literature corresponding to the paramagnetic - 
spin density wave transition in nonsuperconducting LaOFeAs are discussed. In particular, we pay 
attention that upon spin density wave transition there is a relative decrease of the density of states 
on the Fermi level and a pseudogap formation. The values of these quantities are not properly 
described in frames of the density functional theory. The agreement of them with experimental 
estimations becomes more accurate with the use of fixed spin moment procedure when iron spin 
moment is set to experimental value. Strong electron correlations which are not included into the 
present calculation scheme may lead both to the decrease of spin moment and renormalization of 
energy spectrum in the vicinity of the Fermi level for correct description of discussed characteristics. 



Stimulated by the discovery of a new class of high-Tc 
superconductors on the base of LaOFeAs compound [3] , 
a lot of investigations of the electronic and magnetic 
structure of this nonsuperconducting parent compound 
in frames of the Density Functional Theory (DFT) were 
performed [3, 0, I3i S 0- These calculations were 
successful in prediction of not only magnetic instabil- 
ity [H, Q ) but even an exact type of magnetic structure 
of LaOFeAs 0,0]. 

Calculated iron magnetic moment is close to 2 fig 
1^. However, experiments indicate much smaller 
value. Powder neutron diffraction measurements 0] give 
0.36(5) fj,B- Local probe measurements of magnetic 
properties of LaOFeAs such as ^^Fe Mossbauer spec- 
troscopy [lol | together with muon spin relaxation 
indicate the values of ^--^0.35 and 0.25(5) he, corre- 
spondingly. 

The situation when DFT calculations predict larger 
spin magnetic moment in comparison with experimen- 
tal one is rare and known only for few systems (e.g. 
MnSi, ZrZn2 etc.). The inconsistency between experi- 
mental and calculated magnetic moments in these ma- 
terials may be ascribed to spin fluctuations which lead 
to the suppression of magnetic moment ^12] . Neverthe- 
less, LaOFeAs is an outstanding compound even among 
these systems because the ratio locale/ ^-exp is extraordi- 
nary large, ~6, and more importantly since Fe ions in 
simple atomic picture is expected to have S = 2 which 
cannot be easily suppressed by any quantum fluctuations. 

The other known experimental parameters which can 
be compared with their theoretical values are specific 
heat coefficient 7 related to the density of states (DOS) 
on the Fermi level N{Ef) as 7 = ^k^NiEp) and Pauh 
susceptibility x = Mb^(-E'f)- Specific heat coefficient 
can be extracted from the low-temperature behavior of 
the heat capacity. Unfortunately, this parameter is ill 
defined experimentally, i.e. strongly depends on the tem- 
perature range used in the fitting procedure. It was es- 
timated by different groups to be 3.7 mJ/(mol-K2) 0], 
0.9 mJ/(mol-K2) [13] and 0.69 mJ/(mol-K2) Q. How- 
ever, 7 obtained in nonmagnetic DFT calculations 0, 0] 



overestimates the largest experimental value almost in 
two times. The results of magnetic calculations for the 
real striped antiferromagnetic structure [s^ improve the 
situation. They are close to the intermediate experimen- 
tal 7 value. However, one could consider this coincidence 
as accidental since the electronic structure of this antifer- 
romagnetic solution corresponds to the large iron mag- 
netic moment. Susceptibility calculated in frames of non- 
magnetic DFT is 8.5x10"^ emu/mol [1]. At the same 
time the flat region of experimental susceptibility curve 
has the value of ^50x10^'^ emu/mol [l|, 15] - 6 times 
larger than calculated one. 

There is an experimental indication of the formation 
of partial energy gap (or pseudogap) around the Fermi 
level which removes parts of the DOS or few bands from 
the Fermi energy at the phase transition from param- 
agnetic to spin density wave (SDW) state. Direct ex- 
perimental estimations of pseudogap Ep„ is based on 
the results of reflectivity measurements [7|. The pseu- 
dogap where the decrease of optical absorption spec- 
tra is observed at different temperatures corresponding 
to paramagnetic and SDW states is in the range of 
150-350 cm"^ (19-43 meV). The value of the pseu- 
dogap due to SDW formation evaluated from thermal 
transport experiments equals to 210 K [13] which cor- 
responds to 18 meV. Results of temperature-dependent 
angle-integrated laser photoemission study for fluorine 
doped compound indicate ~100 meV pseudogap Il6l . 
whereas high-resolution photoemission spectroscopy [l7| 
gives the pseudogap value above Tc of 15-20 meV. 

The evaluation of the relative reduction of the DOS in 
the vicinity of the Fermi energy due to formation of such 
pseudogap is more indirect procedure. Based on the re- 
sults of susceptibility measurements 15[ one can deduce 
from the values of Pauli-like susceptibility curves at high 
and low temperatures that the change in the N{Ep) does 
not exceed 20%. The same value of —20% one could ex- 
tract from the conductivity curves of Ref. U 

keeping in 

mind that conductivity is proportional to N'^(E) in the 
proximity oi Ep. Analysis of the specific heat at vari- 
ous temperatures [isj suggests 70% reduction of N{E) 
around Ep at paramagnetic - SDW transition. 
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FIG. 1: Band structure of LaOFeAs obtained from nonmag- 
netic calculation in BZ corresponding to enlarged {\^a x 
\/2a X c) unit cell. Grey insert: band structure along S — F 
line for conventional {ax axe) unit cell. Fermi energy is zero. 
For the bands crossing the Fermi level orbital projections are 
marked. 



FIG. 2: Band structure of LaOFeAs for striped antiferromag- 
netic state. Pseudogap energy region is shown by a grey 
shadow stripe. Fermi energy is zero. For the bands cross- 
ing the Fermi level orbital projections are marked. Super- 
scripts mean majority (maj) and minority (min) spin pro- 
jections. There are contributions from both majority and 
minority states to the band dxy. 



Thus for LaOFeAs, one can assume that there is a de- 
crease in DOS near the Fermi energy of tens percent with 
paramagnetic - SDW transition which corresponds to 
pseudogap formation of 20^100 meV. SDW state is char- 
acterized by iron spin magnetic moment of 0.2^-^0.4 /i^, 
specific heat coefficient of 1~4 mJ/(mol-K^) and suscep- 
tibihty coefficient of '--'50x10^^ emu/moL Conventional 
DFT calculations (both nonmagnetic and antiferromag- 
netic) fail in the correct description of these quantities. 

In the present paper we show that fixed spin moment 
DFT calculations with the magnetic moment fixed at ex- 
perimental value can significantly improve an agreement 
with the experiments with regard to the values of spe- 
cific heat, pseudogap and relative decrease of the N{Ef) 
value. 

The calculations were performed within the frame- 
work of the Tight-Binding Linear Muffin-Tin Orbitals 
(TB LMTO) method [11] using the Generafized Gra- 
dient Approximation (GGA), where exchange potential 
was taken in Perdew-Wang form [l^ . Experimentally de- 
termined structure parameters and atomic positions 
for tetragonal phase and collinear striped antiferromag- 
netic order of Fe ions in layer were used. We assumed 
ferromagnetic interlayer interaction due to negligible in- 
fluence of the antiparallel alignment of spins between 
different FeAs layers and for the simplicity of discus- 
sion. The La(6s,6p,5d,4/), Fe(4s,4p,3d), 0(3s,2p,3d) and 
As(4s,4p,4d) orbitals were included into the basis set. 
The integration in the course of the self-consistency iter- 
ations was performed over a mesh of 18 x 18 x 12 k-points 
in the irreducible part of Brillouin zone. We checked that 
such amount of k-points is enough for the precise calcu- 
lation of the Fermi level position Ep and value of the 
density of states at the Fermi level N{Ef)- Fine mesh 
is important due to the nesting bands near Ep. Calcu- 



lations were performed in \/2a x \/2a x c (four formula 
units) unit cell appropriate for description of striped anti- 
ferromagnetic state. Crystallographic x and y axes were 
directed from iron to its nearest iron neighbors and fer- 
romagnetic chains were running along x direction. 

Fig. [1] demonstrates the results of nonmagnetic calcu- 
lations. The band structure agrees with that obtained 
previously 0]. The Fe bands mainly of t2g origin cross 
the Fermi level. One may notice two-dimensional charac- 
ter of the band structure and clear signs of Fermi surface 
nesting in F — X and F — F directions. 

The nesting effect is usually illustrated in the figure of 
Fermi surface. In order to reveal nesting in the simple 
band structure graph one may plot it along 5 — F line 
for conventional a x a x c unit cell; the result is given 
in the insert of Fig. [T] This 5 — F direction corresponds 
to y — r (or X — T) fine for enlarged V^a x \/2a x c 
unit cell. When the unit cell is doubled, left half of grey 
region folds to the right half. The cross of the folded 
bands happens just on the Fermi level. 

The specific heat and susceptibility coefficients 
recalculated from N{Ep) obtained in nonmagnetic 
DFT approach are jnm — 5.3 mJ/(mol-K^) and 
XAfM=7.2x 10~^ emu/mol. That agrees with values cal- 
culated before. 

The total energy difference between nonmagnetic and 
antiferromagnetic (which is energetically more favored) 
states is 116 meV/ (atom Fe) which is in a good agreement 
with the result of Ref. H. Substantial energy gap between 
different magnetic solutions together with the large mag- 
netic moment (see below) shifts the system away from the 
quantum critical point, where spin fluctuations may play 
an important role. Band structure for striped antiferro- 
magnetic state is shown in Fig. O It differs essentially 
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FIG. 3: Band structure of LaOFeAs for striped antiferromag- FIG. 4: Band structure of LaOFeAs for striped antiferromag- 

netic state with the fixed spin moment value 0.36 hb- See netic state with the fixed spin moment value 0.25 hb- See 

also Fig. [2] caption. also Fig. [2] caption. 



from the nonmagnetic picture. In particular, in the vicin- 
ity of the Fermi level inT — X direction there is only one 
band of d^^ character. The other 4 bands are moved away 
from the Fermi level due to the Stoner splitting. In F — F 
direction three bands remain. Two of them [d^y and dyz) 
have the same origin like as in nonmagnetic state. And 
third band originated from c?^2_y2 and d}^^2_^2 orbitals 
which were completely occupied in nonmagnetic case has 
been appeared around Ep- 

One can define the pseudogap as an energy region 
around the Fermi level where number of bands in Fig. [2] 
is essentially decreased in comparison with that in Fig. [T] 
It is natural to define it between maximum of parabola at 
X k-point and minimum of high-lying parabola in F — X 
direction (see grey stripe in Fig. [2]). The pseudogap de- 
fined in such a way is estimated to be 380 meV which is 
much larger than experimental expectations. 

Calculated iron magnetic moment equals to 1.77 /is 
and specific heat coefficient jmag — 0.99 mJ/ (mol-K^). 
This agrees with the experimental estimations of 7 [l3| . 
However, going from nonmagnetic to antiferromagnetic 
phase the calculated value of N{Ef) changes by factor 
of 6 (600% instead of tens percent). 

This inconsistency in the values of magnetic moment, 
width of the pseudogap and too drastic change in N{Ef) 
going from paramagnetic to magnetic state demands the 
explanation. Below by the use of fixed spin moment pro- 
cedure we simulate experimental value of the magnetic 
moment and investigate how N(Ep) and the width of 
the pseudogap is changing upon decrease of the spin mo- 
ment value. 

Band structure corresponding to the fixed spin mo- 
ment value 0.36 fiB is presented in Fig. [3] Now it looks 
very similar to nonmagnetic picture. The remarkable dif- 
ference occurs along T — X direction: d^y and dzx bands 
are spin-splitted and first of them is removed from the 
vicinity of the Fermi level whereas is still crossing 



Ep- There is no dramatic reconstruction of bands along 
Y — T direction. One can conclude that upon transition 
from nonmagnetic to magnetic state the first changes of 
band structure occur along T ~ X line; then the bands 
along Y — T are involving to the formation of magnetic 
moment. 

Fixed spin moment calculation results in a significant 
growth of the specific heat coefficient, which equals to 
IFSM = 2.0 mJ/(mol-K^) in this case in a good agree- 
ment with experiment. The pseudogap which now could 
be defined as it's indicated by grey stripe in Fig. [3] de- 
creases down to 180 meV, which is still larger than ex- 
perimental one. 

It is interesting to note that further decrease of the 
magnetic moment in the fixed spin moment calculation 
(Fig. m leads to even better agreement with experiment 
with respect to the value of the pseudogap. For one of 
the reported values of /i = 0.25 //_b it decreases down to 
130 meV. The specific heat parameter is calculated to be 
2.4 mJ/(mol-K^), and then the change in N{Ef) in non- 
magnetic and magnetic states is only 55% in reasonable 
agreement with experimental estimations. 

Semi-empirical fixed spin moment approach demon- 
strates that the correspondence of experimentally known 
parameters of electronic structure of LaOFeAs is essen- 
tially improved in comparison with conventional mag- 
netic calculations, if the magnetic moment is kept at low 
value ~0.3 fis found in the experiment. However, even 
with reduced spin moment value the calculated pseudo- 
gap (130 meV) remains larger than the experimental one 
(20-40 meV). 

In our opinion, an account for dynamical correlations 
which no doubts exist for d shell of iron may lead both to 
essential reduction of spin moment value and renormal- 
ization of energy spectrum in the vicinity of the Fermi 
level for decreasing of pseudogap. 
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